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Abstract 

The surface free energy is the difference between the free energies for a 
system with open boundary conditions and the same system with periodic 
boundary conditions. We use the quantum transfer matrix formalism to ex- 
press the surface free energy in the thermodynamic limit of systems with 
integrable boundary conditions as a matrix element of certain projection 
operators. Specializing to the XXZ spin 1/2 chain we introduce a novel 
'finite temperature boundary operator' which characterizes the thermody- 
namical properties of surfaces related to integrable boundary conditions. 
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1 Introduction 

The embedding of non-magnetic impurities into spin chains has profound effects on the 
low-energy properties of such systems. Quantities such as the magnetic susceptibility 
will acquire corrections of order 0(L°) which are directly accessible to experimental 
observation (L being the length of the chains). The dependence of these surface con- 
tributions on temperature and magnetic field has been a focus of theoretical work for 
some time now. Apart from field theoretical approaches [7, 8] the exact solution of the 
XXZ spin- 1/2 chain with open boundary conditions [1, 17] has been used to address 
this problem [4, 6, 22]. At zero temperature, the calculation of surface contributions 
to the ground state energy and the magnetic susceptibilities by means of the Bethe 
Ansatz is well established [2, 3, 6, 12, 22]. For finite T, bulk thermal properties are 
known from the thermodynamic Bethe Ansatz (TBA) based on the so-called string hy- 
pothesis [20] . This approach has been extended to compute surface corrections to these 
quantities [4,6,22]. The temperature dependence of the boundary susceptibility ob- 
tained this way, however, disagrees with the field theoretical results [7,8] which, when 
combined with the leading T = contribution from the Bethe Ansatz, have been shown 
to provide the correct low temperature asymptotics [3]. In addition, this approach fails 
to reproduce the correct high-temperature behaviour. Since the Bethe Ansatz yields 
reliable results for the open chains at T = 0, the origin of these problems is likely to be 
the improper treatment of the combinatoric entropy which is a central quantity in the 
TBA. First attempts to extend the TBA for the calculation of the surface contributions 
to the free energy exist [21]. It is still an open issue, though, whether this approach 
can resolve the problems mentioned above. 

In this paper we choose a different approach to obtain the surface contribution to 
the free energy from the exact solution by combining Sklyanin's work [17] on the al- 
gebraic construction of integrable open boundary conditions with the quantum transfer 
matrix approach [15, 16, 18, 19] to the thermodynamics of solvable quantum systems 
in one dimension. 

In section 2 we recall the results of Sklyanin's work [17] as far as they are needed 
for our purposes. Section 3 contains a short review of the quantum transfer matrix 
approach to the thermodynamics of solvable quantum lattice models. In spirit and no- 
tation we follow [9], where more details can be found. Section 4 contains our first 
expression for the surface free energy of a system connected to Sklyanin's reflection 
algebra and treatable by the quantum transfer matrix method. The formula is exposed 
in proposition 4 and shows how to obtain the surface free energy as a Trotter limit 
of an expectation value of a product of certain projection operators which are local 
in the space where the quantum transfer matrix acts. The expectation value is taken 
in a certain so-called dominant eigenstate of the quantum transfer matrix. A second 
expression for the surface free energy as an expectation value of a non-local operator 
in the dominant state is derived in section 5 (see equation (66)). Again a Trotter limit 
must be taken, and a novel 'finite temperature boundary operator' appears. From its 
properties, worked out in the remainder of section 5, we shall see that it is a generic 
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object in the context of the quantum transfer matrix formalism applied to open bound- 
ary systems. From section 5 on we restrict ourselves to the example of the XXZ chain 
with local spins j and only comment on the generalization of the results later in the 
concluding section 6. 



2 Integrable boundary conditions 

Our construction of the finite temperature boundary operator is based on a combina- 
tion of the quantum transfer matrix approach to the thermodynamics of 'Yang-Baxter 
integrable' systems, as reviewed in the following section, with Sklyanin's work [17] 
on the construction of integrable boundary conditions. Sklyanin's construction is valid 
for a fairly general class of integrable systems characterized by ani?-matrix of the form 
R{X,/u) = R(X-/j) e End(V<g>V) (V vector space with dim V = d G N) which not only 
satisfies the Yang-Baxter equation 

7?i 2 (X, A /) J R 1 3(X,V) J R 2 3(A/,v)= J R 23 (A/,v)i?i3(X,V) J R 12 (X, A /) (1) 

but several additional conditions. Namely, R(X) is symmetric in the sense that 

PR(X)P = R(X), (2) 

where P is the transposition map on V (g> V (Px ® y = y <g) x). Furthermore, R(X) is 
unitary, 

R(X)R(-X) = p(X) (3) 
for some complex function p(X), and crossing unitary, 

R tl (X)R tl (-X-2r [ ) = p(X) (4) 



for another complex function p(X) and T| a parameter characterizing the 7?-matrix. The 
superscript t\ in (4) denotes the transposition with respect to the first space in the tensor 
product V <g> V . Similarly, ? 2 will denote the transposition with respect to the second 
space. Then R(X) is a symmetric matrix if 

R V2 (X) =R(X) & R h (X)=R t2 (X), (5) 



which is also assumed by Sklyanin. 

Example. Our chief example will be the 7?-matrix 



R(X) 



(\ 0\ 

b(X) c(X) 

c(X) b(X) 

\0 1/ 



(6) 
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with 

shgj sh(ii) 



This is the well-known 6-vertex model solution of the Yang-Baxter equation (1). It 

l 

2- 



generates the Hamiltonian of the XXZ chain with local spins A. satisfies (2)-(5) 



with 

sh z (X + ri) 

Every solution 7? (A,,//) of the Yang-Baxter equation (1) determines the structure of 
a Yang-Baxter algebra with generators Tp(k),a,$= 1 , . . . , d, through the relations 

R 12 MTi(X)T 2 (m) = T 2 (^)T X {X)R X2 {X^) , (9) 

where 7i(A) = T(k) <g>idy, T 2 (X) = idy ®r(A,). The matrix T (k) of generators is called 
the monodromy matrix. The map defined by 

T a (X) = {T- l )\X) (10) 

is an algebra automorphism, and Sklyanin calls it the antipode. He restricts his con- 
siderations to representations of the Yang-Baxter algebra (9) possessing the crossing 
symmetry 

(T a ) a (k)=®(k)T(k-2r\). (11) 

Here $(A,) is a complex function depending on the representation and T| is the same 
parameter as in (4). 

Sklyanin's construction of integrable systems with boundaries is based on the rep- 
resentations of two 'reflection algebras' and which, for given R(k) satisfy- 
ing (l)-(5), are defined by the relations 

= r 2 { -\ l j)R l2 (l+ l j)T l { -\l)R l2 (l- l j), (12a) 

7? 12 (-X + A /)T 1 (+)fl (^i2(-X- A /-2ri)T 2 (+)f2 ( A /) 

= ^ t \^R l2 (-l-^-2^)rl +)t \X)R l2 (-l+^. (12b) 

We shall call T ( - ) and T ( + ) left and right reflection algebras, respectively. T (~) and 
q~(+) are isomorphic and several isomorphisms are explicitly known [17]. Using (12) 
and (2)-(5) Sklyanin proved the basic 

Theorem 1. The functions ( 'transfer matrices ') 

t(l)=tvT^\l)T ( -\X) (13) 

defined in the tensor product t( + ) ®T ( _ ) form a commutative family, i.e. 

[t(X),t(/j)} = for all X,/jeC. (14) 
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Thus, given a pair of representations of T ( + ) and T (~) the matrix t(X) provides a 
generating function of 'quantum integrals of motion'. Further following Sklyanin [17] 
we shall now review a set of basic results on the representation theory of T ^ and 
T that will be needed in the sequel. 

Proposition 1. Let T^ ± \X) two representations ofT^, respectively, in spaces W ± , 
and T^\X) two representations of the Yang-Baxter algebra (9) in W ± such that 
T^\X) satisfies (11). Then 

T (") (X) = rH (X)?H (X)T^-\-X) , (15a) 
<T = T^\X)V ^\X)T^ a {-X) (15b) 

are representations ofT^ in W ± <S> V7 ± . 

Proposition 1 allows one to construct integrable open boundary conditions for 
known models with L-matrix L(X). These are connected with the simplest representa- 
tions, namely c-number matrices, of T W . Given two representations (X) of T ^ 
in C and an L-matrix representation L(X) (acting on 'a local quantum space') of the 
Yang-Baxter algebra (9), define 

T^(X)=L M (X)...L 1 (X), (16) 
T^+\X)=L L (X)...L M+l (X), (17) 
?( ± \X) = K( ± \X). (18) 

Then, by^ proposition 1, tH(^) = T^~\X)i: H(X)7'(-)- 1 (-X) and <r(+) r (?i) = 
r( + ) r (X)T ( + ) r (X)r( + ) a (— X) are representations of T ^ and, according to theorem 1, 
t(X) = trT^(X)T^(X) generates a commutative family of operators. Moreover, the 
following proposition holds. 

Proposition 2. The generating function t(X) = trT^ + \X)T^\X) with tW(1) de- 
fined by (15)-(18) can be written as 

t(X) = trK^ + \X)T(X)K { -\X)T- l (-X) , (19) 

where T(X) = T^ + \X)T^(X) = Ll(X) . . .L\ (X), and is thus independent of the fac- 
torization (16), (17) ofT(X) into T^+\X) and T^\X). 

Let us now consider any fundamental model associated with a regular 7?-matrix 
satisfying R(0) = P (for a pedagogical review of such type of models see e.g. chapter 
12 of [5]). Its L-matrices have matrix elements L& (X) = R^l(X)ej^, j=l,...L, where 

the ea are matrices with a single non-zero entry 1 at the intersection of the ath row and 
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the fhh column, a, |3 = 1 , ... d. The Hamiltonian associated with a fundamental model 
is 

L-l 

H=(trT(0))- l d x trT(l)\ x=Q = £ ff, 7+1 +ff L1 , (20) 

7=1 

^ = a x (Pi?)^(x)| x=0 . (2i) 

The occurrence of the term Hn on the right hand side means that we are dealing with 
periodic boundary conditions. Using t(X) defined in (19) instead of the 'row-to-row 
transfer matrix' try" (A) as a generating function, we obtain the same Hamiltonian but 
with different boundary terms. 

Proposition 3. Consider the generating function t(X) defined in proposition 2 with 
T(k) a monodromy matrix of a fundamental model and (A) c-number representa- 
tions of the right and left reflection algebras (12) satisfying 

tr/?W (A) = 1 , 7?H (0) = id v . (22) 

Then 

t(k)= + 0(l 2 ), (23) 

where 

* = I,Hjj + i + ^!^ +tro^ +) (0)^o (24) 

7=1 

and Hjk is the same as in (21 ). 

Here the two on-site terms on the right hand side of (24) may be interpreted as 
external fields acting at the boundary sites 1 and L of the chain. In particular, it may 
be possible that both of these terms vanish for certain choices of K^\X) (see next 
example below). Then we are dealing with a chain which is cut at the link between 
sites L and 1, and hence has reflecting ends. Note that, unlike Sklyanin, we have 
assumed that trK^ (X) = 1 for all A, e C. This is not much of a restriction since the 
Hamiltonian turns out to be trivial if trK^ (k) vanishes identically, and otherwise the 
trace can be replaced by 1, since the reflection equations (12) are homogeneous. With 
tr7?(+)(A) = 1 the expression (24) for the Hamiltonian looks slightly simpler, and we 
have t(0) = 1 which will become important later, when we proceed to the quantum 
transfer matrix formalism. 

Example. Diagonal solution of the reflection equations for the XXZ chain [17]. Let 

K(l,Q=I 2 + th(l)cth^)o z , (25) 

where I 2 denotes the 2x2 unit matrix and a z = e\ — e\ the diagonal Pauli matrix. 
Then, up to the multiplication with an arbitrary function, 

*H(X) = *(U~)» tf+\X) = iK(\ + i\,Z+) (26) 



7 



are the unique diagonal c-number solutions of (12) with i?-matrix (6). They are nor- 
malized such as to fulfill (22). Inserting (6), (7), (21) and (25), (26) into (24) we obtain 

2sh(T!)^ = £ [o*0* +1 -G^., +ch(Tl)(0^. +1 - 1)] 

7=1 

+ sh(Ti)[cth(^-)a z 1 +cth^ + )a£] -ch(ri). (27) 

This is the open XXZ spin-^ chain with anisotropy parameter A = ch(r|) and longitu- 
dinal local magnetic fields h T = sh(r|) cth(^ =F ) acting at the left and right boundaries. 
Note that the boundary fields vanish for ^ = ijt/2. 

3 Quantum transfer matrix and free energy 

We recall that every solution R(k,fi) = R(k — fi) of the Yang-Baxter equation (1) de- 
fines an integrable vertex model with row-to-row transfer matrix 

tr T (K) = tr R 0L (X, 0) . . . R Ql (X, 0) = tr L L (I) . . . L x (I) . (28) 

If R(k) is regular, then txT(X) generates a local Hamiltonian with periodic boundary 
conditions via (20), (21). Hence, knowing the eigenvalues of trT(k) means to know 
the eigenvalues of the Hamiltonian (20). 

The row-to-row transfer matrix (28) is not the most convenient tool, however, if 
one wishes to study the thermodynamic properties of the Hamiltonian (20), since, for 
the calculation of the partition function, all eigenvalues of H, and hence all eigenvalues 
of trT(k), are needed. Fortunately, another auxiliary vertex model can be constructed 
by means of the same /^-matrix, whose partition function in a certain so-called Trotter 
limit is equal to the partition function of the Hamiltonian (20) and, asymptotically 
for large L, is determined by a single leading eigenvalue of the associated so-called 
quantum transfer matrix [18, 19]. 

For j = 1 , . . . , L and |3 G C we define the monodromy matrices 

7f>) = ■ ..^(-lX)R fl (X, |) . (29) 

By slight abuse of the usual terminology we shall call them the 'quantum monodromy 
matrices'. They define representations of the Yang-Baxter algebra as is obvious from 
the defining relations (9) and their transpose with respect to space 1 (see [9] for details). 
It is easy to see (e.g. [9]) that 

lim trj N TQ™{0) . . . TQ™{0) = , (30) 

which is the statistical operator associated with the Hamiltonian (20) if we identify |3 
as the inverse temperature 1/7. It follows that the partition function of the L-site chain 
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with periodic boundary conditions is 



CO 



Z L = limtr I ..^(trr^( )) L =£A^(0), 



(3D 



where A n (k) is the nth eigenvalue of the quantum transfer matrix trT®™(X) in the 
Trotter limit N — > °°. By definition Ao(A,) is the eigenvalue with largest modulus at 
X = 0. It must be finitely degenerate in a vicinity of zero. For simplicity we shall 
assume it to be non-degenerate. In any case, the free energy per lattice site in the 
thermodynamic limit turns out to be [18, 19] 



For many models, e.g. for the XXZ chain determined by the i?-matrix (6), the lead- 
ing eigenvalue of the quantum transfer matrix can be calculated for every finite Trotter 
number ./V by means of the algebraic Bethe Ansatz. The corresponding solutions of the 
Bethe Ansatz equations determine auxiliary functions which can be shown to uniquely 
solve certain non-linear integral equations [5,15,16] containing the Trotter number N 
as a mere parameter. In these equations the Trotter limit is easily performed analyti- 
cally, and an integral formula involving the auxiliary functions can be derived for the 
free energy per lattice site (32). Furthermore, it was shown in [9] that the eigenvector 
|*Po) corresponding to the leading eigenvalue Ao(X) of the quantum transfer matrix en- 
codes the complete information about the state of thermal equilibrium and hence about 
all thermal correlation functions. This led to the discovery of integral representations 
of correlation functions of the XXZ chain at finite temperatures [9, 10]. We shall refer 
to |*Po) as the 'dominant state' of the quantum transfer matrix. 

4 The surface free energy as a matrix element 

The derivation of the formula (30) for the statistical operator in the previous section 
relies on a simple generalization of the Euler formula for the exponential function. 
Consider a sequence of operators (X N )^ =l on a finite dimensional vector space that 
converges (in some appropriate norm) to a limit X. Then 



We shall call this the Trotter formula. 

The Trotter formula can be combined with (23) in order to obtain the statistical 
operator for the Hamiltonian H of the open boundary system as a limit. Setting 



/=-rinAo(0). 



(32) 




(33) 



X N = N t(-$j)-l 



(34) 
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and using (23) we find that limN^^X^ = — p\# and therefore, using (33), 

«m,«(-g) =.-!". (35) 

Ne2N 

This observation will serve us as a starting point to express the surface free energy as 
a matrix element of a product of local operators. 

For a pair of c-number solutions K^(X) of the reflection equations (12) let us 
define an operator 

U u {X)=K^\X)P t ^K i r\X), k,£=l,...,N. (36) 

Proposition 4. The surface free energy, i.e. the difference between the free energies 
of the open system with Hamiltonian (24) and the corresponding periodically closed 
system with Hamiltonian (20) is given by 

AF= lim -rin(Vo|n E (-g)...n 3 ^(-g)|V >, (37) 

Nem 

where |*Po) is the dominant state of the quantum transfer matrix \xT@™ (X) (see (29)). 
Proof. Let us define 

T- k (X)=R- kL (X,0)...R- kl (X,0), (38a) 
T- k (X)=R l - k (0,X)...R Ll (0,X) (38b) 

for k = 1, . . . ,N (remark: then T^(X) = Tj~ (X) for unitary R matrices with p(X) = 1*). 
It follows that 

t(X) = tr l K{ + \X)T l (X)K ( [ -\X)T l (-X) 
= tr u K[ +) (X)PaTi (X)K{- } (X) T- x ( -X) 
= tr u K[ + \X) {m)Pa) ^(XWrt-X) 

= tr I5 ^ +) (X)4^- ) (X)r|(X)r I (-^) 

= tr u U- 1 - 2 (X)mX)T- 1 (-X). (39) 

Hence, 

= tr L.JV n i2 n 34 • • • n N^lN T h r A^T(l) • • • T { {~Tj) T l (I) 

= tr L ..^ n H n» . . . n^zf r >) . . . r L er » . (40) 

*Henceforth we simply assume that p (A.) = 1 . This can usually be achieved by properly normalizing 
the /^-matrix and is certainly true for the /^-matrix (6) of the XXZ chain (see equation (8)). 
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For the last line compare equation (24) of [9]. For simplicity we left out the argument 
— P//V of the operators n here. Using (35) we obtain the partition function Zl of the 
L-site open system as 



z L = lim tr I ..^n I ,(-|)...n^(-|)(trr^ ( o)) L 



A'e2N 



d N -i 



= ^L I (^|n I2 (-|)...n^(-|)|^ M )A^(o), (4i) 

Nem 11=0 

where (*P„| and |*P„) are the left and right eigenvectors of the quantum transfer matrix. 
Comparing (31) and (41) we conclude that 

fii^ = Wn E (-g)...n^(-g)|v > (42) 

which completes the proof of proposition 4. □ 
Let us examine the structure of the operator n, equation (36). Its components are 

ng(X) = ^ + )(X)«^-)(X)f. (43) 
This suggests to introduce row and column vectors | — ) and (+| with components 

a p\-) = k(~\X)1 (+\^ = K^\X)l (44) 

Then 

ng(X)= p8 |->(+r=(|->®(+|)g. (45) 

As usual we leave out the tensor product sign and write II (A,) = |— )(+|. From the 
definitions of |— ) and (+| it follows that |— )(+| is proportional to a one-dimensional 
projector, 

n 2 (X) = (+|-)n(A), (+|-)=tr*< + >(X)*<->(a,). (46) 

Example. Let us illustrate our result with the example of the diagonal '^-matrices' 
(26) of the XXZ chain. For these matrices 

i- >= s(?4d^) (shr+X) ' ' 0,shR "" X)) ' <47) 

(+l = 2shft + )'h(A + 11 ) < Shft+ + A + ll) ' - - Shft+ -"- 1l) )'' (48) 

Hence, 

n(X) = |-)(+| 



2sh(^+)sh(^-)ch(A)ch(A + ri) 



/sh(^-+A)sh(^+ + A + ri) sh(^--A)sh(^++A + ri)\ 





Vsh($- + X)sh($ + -X-Ti) sh(Zr-X)sh(Z+-X-T\)/ 



(49) 
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and 



(+|_) = 1 + cth(^-) Cth(£+) th(l) th(l + T1) . 



(50) 



As we remarked above, open boundaries with vanishing boundary fields are realized 
for^+ = =i%/2. Then sh(i7l/2 + X) = ich(X) and thus 



n(X) = - 



(\ 
o 



V 1 










1\ 






(51) 



which is a ^-independent true projection operator. 



5 The finite temperature boundary operator 

In proposition 4 we have expressed the surface free energy as an expectation value of a 
product of certain projection operators n in the dominant state of the quantum transfer 
matrix. The projectors n are local in the space where the quantum transfer matrix 
acts. Having in mind that, within the algebraic Bethe Ansatz, the dominant state is 
constructed by acting with certain creation-operator-like monodromy matrix elements 
on a 'Fock vacuum', we would like to express the product of projection operators in 
equation (37) as well in terms of elements of the monodromy matrix. This is possible 
by means of an inversion formula obtained in [11, 13]. 

Instead of working out all steps in full generality we shall stick from now on with 
the example of the XXZ chain. The final result will be of sufficiently general appear- 
ance, and we shall indicate how to move on to other systems. In case of the XXZ chain 
the auxiliary vector space V is two-dimensional, and the quantum monodromy matrix 
(29) becomes a 2 x 2 matrix when expressed in terms of the L-matrices 

*„(M = H +* ( VH *tt C( Vu* 2) , (52) 
V c{X-n)e n \ b(\-n)e„[+e n $J 

7 {„ u _ (e n \+b{n-X)e n l c{^i-X)e n \ \ 

L ^)-{ c{n-X)e n \ " b(M-X)e n \+e H l) > (53) 

namely, 

r^» = 2 Ar (-|,x)x^_ 1 (xJ)...2 2 (-|,x)x 1 (x,|). (54) 

The aforementioned inversion formula [11, 13] applies to the monodromy matrix 
corresponding to the usual row-to-row transfer matrix, not to the quantum monodromy 
matrix. Fortunately however, in case of the XXZ chain the different types of mon- 
odromy matrices are simply related by means of the crossing symmetry. 



(55) 
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where <Jj is the Pauli matrix a y = (9 q) acting on site j. Setting 

S = G y N c y N _ 2 ...G y 2 , a(X) = ^(-$-X), (56) 

we find that 

^^^ = s^(M-fiW-i(x,|)...x2(M-fi)^i(^S)5 

= Sr(^i,...,^v)S, (57) 

where , . . . , ^v) =: T(k) is the inhomogeneous row-to-row transfer matrix with 

inhomogeneities 

T| - I if J even 



(58) 



if 7 odd 

The right hand side of (57) admits the application of the inversion formula [11,13] 

rl-l -, r N 



7=1 j=n+l 



(59) 



The left hand side is the rescaled quantum monodromy matrix. We abbreviate it as 



TQ™(1) ~ otm fA(X) B(X) 



a(X) 



\C(X) D(X) 



(60) 



Taking into account that 

Se n p = (oy) n+l "e n l(oy) n+l l S 2 = \ (61) 

we obtain the following inversion formula for the rescaled quantum monodromy ma- 
trix (60). 



Proposition 5. 

r n-l 



*4= 



N 



ntrf^^) [(oyr + 'fQ™(k n )(&r + % n ktq™^) 

7=1 j=n+l 



(62) 



This has an immediate application to our problem of expressing the projection 
operator IT, equation (36), in terms of monodromy matrix elements, 

n 2 „-l 2 „(-|) = IT trf2™(y x(-|) J] , (63) 

J=l l j=2n+l 

x(X) = trK { -\l)f QTM (-l)K (+ \l)a y f QTMt (X + T])a y . (64) 
We shall call x the 'finite temperature boundary operator'. Equation (63) implies that 

M-DM-!) . . .iw(-g) = x^ (-£ ) . (6 5) 

Inserting the latter equation into (37) we can formulate the following theorem. 
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Theorem 2. The surface free energy of the open XXZ chain is determined by 

AF = i in L- rln — mm) » (66) 

where \ {X}) is the unnormalized dominant state of the reseated quantum transfer ma- 
trix trT@™ \X), equation (60), and x is the finite temperature boundary operator as 
defined in (64). 

Recall that within the scheme of the algebraic Bethe Ansatz the vectors \{k}) and 
({A,} | are given by 

\{X}) = B(Xi) . . .B(Xn)\0) , ({X}\ = (0\C(X» ) . . .C(Xi) , (67) 

where B(X) and C(X) have been defined in (60), where |0) = [(J) <8> (O)]®^ 2 is the 
staggered pseudo vacuum of the quantum transfer matrix with dual (0|, and where 
{Xj} N j=\ is the solution of the Bethe Ansatz equations 



Sh(l; + |) Sh(^ - i +T|) 1 T ^ S h(\j - Xj, + tl) 

-iW>-»*-u' ( ' 



S h(^-|)sh(^+|-Tl) 



j = 1, . . . ,N/2, that maximizes the modulus of the leading eigenvalue of the quantum 
transfer matrix at A = 0. 

Example. As an example let us again consider the XXZ chain with reflecting ends, i.e. 
vanishing boundary fields. We have to set ^ = irc/2 in (25), leading to 2K^(X) = 
K^~\X) = I2, the 2x2 unit matrix. Inserting this into the definition (64) of the finite 
temperature boundary operator we obtain 

x(X) = \trT Q ™{-X)o y T QTMt (X + x\)<5 y (69) 
= ^(A(-X)D(X + r[)-B(-X)C(X + r\)-C(-X)B(X + r\)+D(-X)A(X + r[)). 

In the remaining part of this section we would like to demonstrate that the finite 
temperature boundary operator x(X) is generically connected with the reflection alge- 
bra (12). Our basic result is 

Proposition 6. With any c-number representation K^ + \X) of the right reflection alge- 
bra (12b) the matrix 

<T (+)'(*.) = f Q ™(-X)K {+ \X)o y f QTMt (X + r[)o y (70) 

as well defines a representation of (12b). 
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Proof. We only sketch the proof. The idea is to show that 

T (+ \K):=T QTMt (-X), (71) 

which clearly defines a representation of the Yang-Baxter algebra (9), satisfies (11), 
and then to employ proposition 1 . In order to show that (X) satisfies (11) one may, 
for instance, use (57), (60) and apply to the elementary L-matrix (52) equation (55) 
and the general relation 

for the inverse of a monodromy matrix associated with the 7?-matrix (6). Notice that 
det q T(X) in (72) denotes the quantum determinant (see e.g. [17] or chapter 14 of [5]). 
For (k) the function f>(X) in (1 1) turns out to be 



M-x+|-tDM-X-|) i * 

i,(-X-|+T!)H-X+|)J 



Knowing that T^ + '(X) satisfies (11) we apply proposition 1 to the special case, when 
T = K^ + \X) is a c-number representation of the right reflection algebra (12b). 

Then also K( + )*(X) is a c-number representation of (12b), which can be seen by taking 
the transpose of (12b) and using the symmetry of the 7?-matrix. Taking into account 
that 

r(+) = ygra-i^T!) = y T QTM^ l + 7])a y 

{ ' fl(-X) ^K-X)deV^e™(?l+2r|) , 

the proposition follows. □ 

Proposition 6 has a number of simple but important consequences. 

Corollary 1. For any pair of c-number solutions K^>(X) of the reflection equations 
(12) the following holds. 

(i) The finite temperature boundary operator x(X), equation (64), defines a com- 
muting family of operators, 

[t(X),t(a*)] = foralll^eC. (75) 

( ii) x(X) can be diagonalized by means of the algebraic Bethe Ansatz. 
( Hi) The matrix element 

N/2 

*(vi,...,v,) = (WinCim) (76) 

7=1 

is a symmetric function of the parameters Vj G C, j = 1, . . . ,N/2. 
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6 Conclusions 

We have expressed the surface free energy for the XXZ chain in the thermodynamic 
limit as Trotter limits (37) and (66) of expectation values of certain operators in the 
dominant eigenstate of the quantum transfer matrix. 

In (37) the operator is local and is the product of certain projection operators n 
defined in (36). This formula is not only valid for the XXZ chain with reflecting ends, 
but was derived for a rather large class of open boundary systems, specified by a pair of 
c-number solutions of the reflection equations (12). If the corresponding i?-matrix of 
'difference form' and satisfying (2)-(5) is regular we can always associate a quantum 
transfer matrix with it, and the projection operators (36) are well-defined. 

The operator in (66) is non-local and is a simple quadratic expression in terms of 
the matrix elements of the rescaled quantum monodromy matrix (60). We called it 
the finite temperature boundary operator. In proposition 6 and corollary 1 we worked 
out its relation to the reflection algebra (12). We derived (66) for the case of the XXZ 
chain with boundary conditions specified by any c-number representation of the reflec- 
tion algebra (12). Still, we used only a few rather general properties of monodromy 
matrices connected to the 7?-matrix of the XXZ chain, namely the crossing relation in 
the form (55) and the existence of a quantum determinant, which allowed us to apply 
the inversion formula (72). It seems rather likely that a finite temperature boundary 
operator can be defined for all solvable systems having these additional properties. 

The explicit calculation of the surface free energy from our formulae remains a 
challenging task, since the calculation of matrix elements, such as in (37) and (66), 
within the algebraic Bethe Ansatz is, in general, involved. In the light of the results 
obtained e.g. in [10, 14] there is, however, hope to accomplish this challenging task in 
the future. 
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